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ABSTRACT 


.  The  present  paper  considers  the  stochastic  difference  equation  Y^  « 
where  and  Qn  are  respectively  random  d  *  d  matrices  and  random  d-vectors,  and  ob¬ 
tains  some  reasonable  sufficient  conditions  on  and  Qn  under  which  Y£  converges 
in  distribution.  In  addition,  a  particular  model  is  examined  when  d  *  2.  in  which 
the  asymptotic  independence  of  #  and  results  in  a  characterization  of  the 

Gamma  distribution. 


1.  INTRODUCTION 


In  this  paper  we  study  the  limit  distribution  of  the  solution  Yn  of  the 
difference  equation 

(1.1)  VV-I 

where  M_  are  random  d  *  d  matrices,  and  0  and  Y_  are  d- vectors  (where  d  could  be 
considered  infinite,  unless  otherwise  stated) .  Throughout  we  take  the  sequence  of 
pairs  (Mn,  Q^,  n  2  1,  to  be  independent  and  identically  distributed.  Equation  (1.1) 
first  came  to  our  attention  in  a  paper  by  Bernard,  Shenton,  and  Uppuluri  [l],  in 
which  it  was  used  as  a  model  for  the  distribution  of  radioactive  material  in  the 
bone  structure  of  humans.  Since  then,  we  have  seen  it  arise  in  a  variety  of  other 
contexts  (see  Soloman  [9],  and  Cavalli-Sforza  and  Feldman  [2],  for  examples) . 

The  asymptotic  behavior  of  (1.1)  is  examined  in  [10]  by  Vervaat  in  the  special 
case  where  d  »  1.  A  variety  of  conditions  are  given  for  the  convergence  in  distri¬ 
bution  of  (1.1).  In  [5] ,  Kesten  establishes  a  reasonably  general  condition  under 
which  (1.1)  converges  for  the  cases  where  d  >  1.  In  [5],  it  is  shown  that  if 


(1.2)  E(ln|  |m|  !)*««• 


then  there  is  a  constant  a  where 

(1.3)  a  s  lim  1/n  ||M,  ...  M  |J  a.s., 

n-*» 

where  for  ad«d  matrix  M  and  a  d-vector  x,  we  define 

I 

(1.4)  | |m| {  *  max  |Mx|,  and 

|x|  -1 

d  u 

(1.5)  |x|  -  (|  xfp. 

i-1  , 

We  then  have  in  [5]  that  if 

-  1  - 


I 


(1.6) 


<  », 


i)  a  <  0  in  (1.3),  and 
ii)  there  exists  a  B  >  0  where  EjQ^1 


then  (1.1)  converges  in  distribution. 

The  first  objective  of  this  paper  is  to  establish  criteria  other  than  (1.6) 
for  the  convergence  of  (1.1).  While  the  conditions  established  will  show  that 
(1.6ii)  can  be  weakened,  that  is  not  the  aim  of  this  paper.  Instead,  the  criteria 
attempt  to  bypass  what  we  perceive  to  be  the  major  difficulty  for  applicability 
of  [5],  that  is,  the  extreme  difficulty  in  the  determination  of  o  in  (1.3). 

In  [7]  and  (10],  the  respective  authors  investigate  possible  limiting  dis¬ 
tributions  of  (1.1)  in  the  case  when  d  ■  1,  and  achieve  some  partial  results  in 
classifying  the  limiting  distributions.  The  second  objective  of  this  paper  is  to 
examine  what  criteria  on  the  model  (1.1)  for  d  >  1  will  yield  a  limit  behavior  in 
which  the  various  components  of  the  vector  Yr  are  asymptotically  independent  for 
a  special  model  in  which  only  a  one-way  reapportioning  of  material  is  allowed.  We 
find  that  in  most  cases  that  non-trivial  asymptotic  independence  is  in  general 
achieved  only  if  the  limit  distribution  of  the  components  is  a  gamma  distribution. 
As  such,  we  arrive  at  a  slight  generalization  of  (Lukacs  [6]) . 


2.  CONVERGENCE  OF  THE  DIFFERENCE  EQUATION. 


From  (1.1),  it  is  easy  to  see  by  iteration  that 


(2.1) 


Yn  * Jj  Mn  * ' '  Mi*lQi  *  Mn 


which,  for  given  Yn  has  the  same  distribution  as 


2.2 


I 


He  shall  establish  that  in  the  case  when  ...  ■*  0  a.s.  exponentially  fast,  under 

weak  conditions  on  Yg  and  (2.2)  will  converge  almost  surely,  independent  of  the 
distribution  of  Yg.  Then  a  variety  of  conditions  which  should  be  reasonably  easy  to 
check  in  particular  cases  will  be  given  which  are  sufficient  to  guarantee  that 
Mj  ...  Mn  converge  to  zero  exponentially  fast. 

For  convenience,  for  any  matrix  M  and  vector  Y,  we  will  use  the  typical  notation 


(2.3)  (m^)  =  M  ,  (yA)  2  Y. 


Also,  let  us  use  the  following  notation,  given  matrices  Mj,  ...  Mr,  and  vector  x; 

(2.4)  M(n)  »  1^  ...  Mn, 

d 

(2.5)  1*1.  •  I  l*±l • 

i=l  1 

(2.6)  I|m||.  1  I^L* 

We  now  establish  the  main  theorem  in  this  section,  which  describes  sufficient 
conditions  for  the  almost  sure  convergence  of  (2.2). 


THEOREM  2.1:  If  there  is  a  0  <  X  <  1  where  either  i)  X"n  sup  (T  lm.  .(n)|)  ■*>  0. 

lsjsd  i=l  x'3 

a.s.,  E(ln|Q_  |  )*  <  »,  and  |Y_1  <  •,  a.s.,  or  ii)x’nsup  (f|m  /n3|)  0,  a.s., 

isisd  j-1  1,3 

E(ln(sup|q,  ,  |)+)  <  »,  and  sup|y„  . |  <  •,  a.s.,  then  R  converges  almost  surely, 
isisd  1,1  lsisd  0,1  n 


PROOF:  The  proofs  that  condition  i)  and  condition  ii)  are  sufficient  for  the  almost 
sure  convergence  of  Rr  are  similar,  and  as  such,  it  will  only  be  established  for 
condition  i). 


-  4  - 


We  will  establish  that  is  a  Cauchy  sequence  with  respect  to  the  infinity 

non,  almost  surely.  That  is,  that  for  all  e  >  0,  there  is  an  N  where 


(2.7) 


P(  lRn"RJ«  *  €*  n  >  N)  >  i_e- 


will  ignore  the  term  ...  M^Yq,  since  it  is  clearly  negligible.  Then  assuming 


n  >  m. 


(2.8) 


Ivj-  “JJ'n.i  -  '..i1 


!  i  t  l-1^*1,IUJ  kl 

sffl+1  i-1  k-1  x»* 
j-m+1  k-1  1-1  Xj*  >’K 


(2.9) 


C.  «  (ui:  n  sup  (X"n|  li  *  «># 

N  n-N  lsfcsd  i-1  4,K 


we  have  from  (2.8)  that 


n  d 


(2.19)  \K-\L  s  !C  C«)  I  U V ’  U«  vl 

n  M  cm  )-*♦!  k-1  J’k 


w  f  I  c  S  \\V%  kl- 

Sn  i-1  k-1  i-1  1,K  3,K 


where  IA(x)  is  the  indicator  function  of  a  set  A. 
Thus,  from  (2.10),  we  have  that 


(2,11) 


P(lRn-Rj„  >  e,  some  n,m  >  N) 

*  P(cN)  ♦  P(  I  \  vl  >  O 

N  j-N+l  k-1  3* 

«  Pfo)  ♦  PC  I  eA,"1|Qj.  >  C). 

N  j-N+1  3 


By  leana  2.2,  to  be  presented  later,  we  have  that  t  X  | Q - J  converges  almost 

i-1  1 

surely  if  and  only  if  EOnjQjJ^)*  <  •.  Prom  condition  i)  we  have  that  there  exists 
a  X,  0<X<1,  such  that  we  have  that  for  sufficiently  large  N,  P(CR)  >  l-e/2  for  all 
n  >  N.  Thus,  it  is  easy  to  see  from  (2.11)  that  for  sufficiently  large  N, 


(2.12) 


pcIV 


>  e,  some  n,  m  >  N)  <  c/2  ♦  e/2  *  e. 


which  establishes  the  theorem.  □ 

In  [10],  Verwaat  establishes  the  lemma  cited  above  (which  is  actually  more 
general  than  the  current  requirements),  which  we  now  state  for  the  sake  of  complete 
ness. 

LEMMA  2.2;  For  <Xn> ,  {Yn>,  i.i.d.  random  variables,  where  -■»  <E(ln|X1 1)<  0,  then 
n  i- 1  ^ 

7  (n  Xj)Y,  converges  a.s.  if  and  only  if  E(ln|Y, |)  <  •. 

i-1  j-1  3  1  1 

We  will  now  establish  several  criteria  which  are  sufficient  for  the  exponen¬ 
tially  fast  convergence  of  MjMj  ...  M^.  To  do  this,  we  will  appeal  to  a  general 
lemma  about  functions  on  the  matrices. 

LEMMA  2.3:  I£  (An)  is  a  sequence  of  i.i.d.  random  matrice,  and  f  is  some  function 

such  that  f(AjAj)  S  f(Aj)  f(A^)  and  E(ln  £(AX))<  0,  then  there  is  a  X,  0<X<1  where 

P(f(  D  A.)  >  Xn  i.o.)  -  0. 
i-1  1 

PROOP:  By  assumption,  we  have  that 


(2.13) 


n  n  n 

In  f(  n  A.)  &  ln(  H  f(A.))  •  7  In  f(A.), 
i-1  1  i-1  1  i-1 


Thus,  if  we  define  6  ■  -E  In  f(Aj),  we  have  for  any  y  where  0<y<6  that 


(2.14)  P(  ?  In  f(A.)  >  -yn  i.o.)  «  0. 

x  * 


As  such,  from  (2.13)  we  have  that 


(2. IS)  P(f(  H  A.)  >  e"'Tx.o.)  £  P(  I  In  f(A.)  >  -vn  i.o.)  -  0. 

i-1  1  i-1  1 

We  get  that  for  all  X  >  e "  the  assertion  is  true.  □ 

Now  we  establish  the  first  sufficient  criteria  for  the  convergence  of 
which  should  be  reasonably  easy  to  verify  in  certain  cases. 

THEOREM  2.4;  If  E  to  (UmJIJ  <  0,  E  InClQjlj*  <  and  l*0L  <  “  2*1  R„ 

converges  almost  surely. 


PROOF :  By  conditional)  of  Theorem  2.1,  we  need  only  establish  that 
X”n  sup  (  l  |m,l?J||  ♦  0,  a.s.  This  follows  quickly  from  Leraaa  2.2  by  noting 


l 

lSjSd  li«l  *3 
that  for  all  x  where  |Bx|m  *  0, 


(2.16) 


|ABx| 


lABxj, 

|bx|4 


IbxI.  s  ||a||Jbx|„, 


where  the  inequality  follows  from  the  definition  of  theorem  (2.6).  Thus  we  have  that 
Hab|L  *  ||a|  I,  ||B||  .  As  such,  from  the  assumptions  and  lemma  2.2  we  have  that 


(2.17) 


H\\  n  m. j |  >  X11  i.o.)  -  o. 
i-1  x 


This  yields  that  for  X  <  X^  <  1,  we  have  that 


(2.18)  x:n| |  n  M.|l  ♦  0,  a.s.  . 

1  i-1  1 

As  can  be  easily  established  from  (2.6),  we  have  that 

(2.19)  II  H  M.JI.  •  sup  (  f  |m  (”5|), 

i-1  1  lijsd  i-1  1,3 


and  the  theorem  is  completed. 


While  the  condition  E  fcn  1 1 j  |  ^  <  0  is  not  as  formidable  a  condition  to  check 
as  at  first  it  appears,  since  formula  (2.19)  gives  some  hope  of  verifying 
E  tnllMjU^  <  0,  nonetheless  it  is  perhaps  too  difficult  for  certain  models.  The 
following  theorem  gives  a  simpler  condition  to  check. 

Using  the  matrix  notation  e|m]  »  (E|m,  . j) ,  we  have  the  following  theorem. 

THEOREM  2.S:  If  d  <  -JIeImJ  <  1,  E  inClQ^J*  <  «,  and  \yq\„  <  -,  a.s., 
then  Rn  converges  almost  surely. 


PROOF:  Once  again,  by  condition  i)  of  theorem  2.1,  we  need  only  establish  that  for 

0  <  X  <  1,  X"n  sup  (  f  -*■  0,  a.s.  Since  ||e|m^|}|  <  1,  we  can  choose  a 

Ujsd  i-1  1,3 

X  and  X^  where  ||e|M^|  |  <  Xj  <  X  <  1.  From  (2.16),  we  have  that 

(2.20)  II  n  e|m.|  II.  s  H  |  |e|m.|  \\m  <  x  n. 

k«l  K  k=l  *  1 

Also,  we  have  from  the  independence  of  the  matrices 


(2.21)  '  «IJ,  I  Sk!1E**l.i.1t,EK.k,jl 


which  alongwith  (2.19)  and  (2.20)  yields 

(2.22)  llEjMj  ...  Mn|  II.  S  II  “  *1^1  II.  <  Xj®. 


Also,  we  have  for  X  >  Xj  that 


P(  ux-n  sup 
n-m+1  lSjSd 


C  l  kC?|)  s  l  P(X-nsup  (  flm.^l)  >  e) . 
i— 1  1,3  n-m+1  isksd  i-1  1,3 


(n)i 


s 


I*'VlB  sup  (  f  | 
n-m+1  lsjsd  i-1 


i  <n) 


l). 


(2.23) 


E  sup  l  |m.^|  *  l  E  £  |m.^|  *  d  sup  l  Elm,^ 

Ujsd  |i»l  1,3  j-1  i-1  1,3  lSjSd  [i-1  1,3  j 


d  llEjMj...  Mn|  || ., 


we  get  from  (2.23)  that 


(2.25) 


u  X‘n  sup  I  4  I  >£  X"1e“1d||E|M1  .. 

n-m*l  lsjsd  i*l  *3  n=m+l  1 


-lXln 


s  I  d  e 
n-m*l 


where  the  second  inequality  follows  from  (2.22). 
As  such,  for  m  sufficiently  large,  we  have 


(2.26) 


•  d  ^ 

P  u  x'n  sup  l  lm.^1  >  e  <  e, 

n«m+l  l£j£d[i-l  ,3 


which  conqpletes  the  proof.  □ 

We  now  introduce  another  functional  on  matrices  which  is  very  similar  to  the 
spectral  radius,  first  introduced  by  Dobrushin  [3]  in  the  case  of  stochastic 
matrices.  For  any  d*d  matrix  P  ■  (Pjj)»  let  us  define  an  auxiliary  matrix 
P  -  (PAj)  by 


(2.27) 


1  if  i-1,  j-1 

0  if  i  >  l,  j-1 

d 

1-  l  p. .  if  i-1,  j  >  1 
i-1  13 

Pi-1 j-1  if  i  >  l»  *  >  1 


Then  we  define  our  functional  6(P)  by 


(2.28) 


*VP  1  fPji”  Pjkl* 

,k£d*ll j-1  l  31  3KJ 


LBMA2.6:  For  any  d*d  matrices  P,  Q,  where 


’  d  f  d 

*«P  l  IPhI  <  ®.  sup  l  |q..|  <  -,  6(PQ)  <  5(P)  5(Q) 

WjsdU-l  3  isjsdU-i  3 


PROOF:  The  functional  6  is  very  similar  to  that  of  Dobrushin  [3] .  An  examination 
of  the  proof  that  for  stochastic  matrices  P',  Q',  we  have  «(P*Q’)  *  «(P')6(Q’)» 
(see,  Isaacson  and  Madsen  [4]),  yields  that  the  crucial  properties  are 


(2.29) 


d*l 

i)  l  p..  -  1.  for  all  i 
j»l 


d+1  d+1 

111  L  lV  -  v  -jli  w 


all  i,  k. 


Since  both  are  satisfied  by  the  definition  of  Pt  Q,  the  proof  is  complete.  □ 

It  should  be  noted  that  6(*)  and  J|*||„  are  very  similar.  In  fact,  if  P  is  a 

non-negative  matrix,  then  it  is  easy  to  show  that  6(P)  *  | |p| 1^.  However,  if  p  is 

allowed  to  have  negative  values,  it  is  possible  that  6(P)  <  J  Jp| or  1 1 p| <  $(p). 

Because  of  this  similarity  in  behavior,  it  is  easy  to  see  that  a  theorem  similar  to 

theorem  2.4  can  be  established  for  fi(ML)  by  methods  identical  to  those  used  in 

d 

theorem  2.4,  (we  need  only  notice  that  from  (2.28)  we  have  6(P)  e  sup  (  J  p..  ♦), 
d  d  lSisd  j«l  J1 

4(P)  i  sup  (  l  p  "),  which  yields  that  26(P)  fc  sup  (  £  |p..|)).  As  such,  we 
ISiSd  j-1  J  lsisd  j=i 

state  without  proof  the  following  theorem. 


XjgggJLL.  ilE  ln6(Mj)  <  0,  E  lndQ^,)*  <  «,  and  |Y0|w  <  -  a.s. ,  then 
converges  almost  surely. 

It  should  be  noted  that  theorems  2.4,  2.5,  and  2.7  were  applications  of  con¬ 
dition  i)  of  theorem  2.1.  It  can  be  easily  seen  that  by  considering  M^',  the  trans¬ 
pose  of  Mj,  that  theorems  similar  to  2.4,  2.5  and  2.7  can  be  established  using  con¬ 
dition  ii)  and  virtually  identical  proofs  to  the  theorems  already  established  in¬ 
stead.  However,  both  the  statements  and  proofs  of  these  theorems  will  be  omitted  in 
this  paper. 


-  10  - 


GB  FO 
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In  [7]  and  [10],  the  respective  authors  investigate  possible  limiting  dis¬ 
tributions  of  (1.1)  for  the  case  of  d*l,  and  achieve  some  partial  results  in  clas¬ 
sifying  possible  limiting  distributions.  For  d®2,  we  will  examine  for  a  special  mo¬ 
del  the  conditions  under  which  the  two  components  are  asymptotically  independent. 

The  model  to  be  considered  is  the  difference  equation 

where  (V  *s  101  i.i.d.  sequence,  and  V _,wn»Qn  ar®  independent  of  each  other 

for  all  n.  This  represents  a  one-way  flow  storage  model,  in  which  at  step  n  new 
material  is  added  to  component  one  via  Q^,  material  is  transferred  from  component 
one  to  component  two  via  l~Vn>  and  material  is  lost  from  the  system  from  component 
two  via  1-W  • 

fl 

For  the  model  given  in  (3.1),  if  we  let  ^  X*  ***<*  let 


(3.2) 


isY.  ♦  ttY- 
<Ks.t)  *  E(e  A  z) 


4»(»)  ■  Ee  , 


then  it  is  easy  to  verify  from  (3.1)  that 

(3.3)  <|>(s,t)  «  *(fl)B*(sV1  ♦  td-Vj),  tWj) . 

Since  Yj  and  Y2  are  independent  if  and  only  if  <Ks,t)  ®  .  w©  have 

from  (3.3)  that 

(3.4)  *(s.t)  -  <>(s)  E^(sV1  ♦  td-V^.O)  EHO.tWj). 

Also  from  (3.3)  it  can  also  be  sl^wn  that 


(3.S)  *(s,t)  •  tfs)  EtfsVj.O)  E^tCl-V^  .tWj) 


Ks)  B+CsV^O)  E^tCl-Vj)  ,0)  E^O.tWj). 


If  we  let 


(3.6)  A  ■  (s:  4»(s)  «  0  or  Et(0,sWj)  e  o), 

and  if  Ac  is  dense,  then  by  equating  (3.4)  and  (3.5)  one  can  see  that  for  all 
(s,t)  €  AcxAc  that 


(3.7)  EKsVj  ♦  tU-Vp.O)  »  E^sVj.O)  E^td-Vj)  ,0) . 


Then,  by  continuity  of  characteristic  functions,  we  get  that  (3.7)  holds  for  all 
(s,t)  cR  *  R.  Thus,  under  conditions  sufficient  to  guarantee  that  Ac  is  dense,  we 
get  from  (3.7)  that  Yj,  Y2  are  independent  if  and  only  if  for  Vj  independent  of 
Yl*  V1Y1  (l-V^JYi  are  independent. 


In  the  following,  we  will  say 


(3.8)  X  ~  r(X,g)if  P(X  s  x)  - 


0  if  x  S  0 


/„  t(«‘1X(Xy)9‘1e"Xydy  if  x  >  0, 


(3.9)  X~B(a,8)if  P(X  S  x)  - 


0  if  x  s  0  or  x  *  1 


/x  B(a,0) ”1  ya_1(l-y)B'1dy  if  0<x<l. 
0 


The  main  result  of  the  section  is  the  following  theorem. 


THEOREM  3.1:  For  model  (3.1),  if  V.  and  Y,  are  independent  then  V.Y.  and 


(l-VjjYj  are  independent  if  and  only  if  one  of  the  following  six  conditions  are 


Also,  it  is  dear  P(U  -  W  <  bj  -  b2  ♦  e)  >  0,  for  all  e  >  0. 

Since  U(U  -  IQ  U  -  W  are  independent,  we  have  that 

(3.15)  P(U(U  -  W)"1  >  b3(bj-b2+e)_1)  -  P(U(U  -  W)'1  >  b3(b3-b2+e)-1|U  -  W  <  b 
And  U  -  If  <  b3-b2-fe  implies 

(3.14)  (J  <  b3+e  and  W  >  b2-e. 

As  such,  we  have  that  U  -  If  <  b3~b2*c  implies 
CJ.1S)  0(U  -  w)’1  >  »  .  . 

i  N 

Thus  P(U(U  -  W)  >  4£  |U  -  If  <  bj-b2+e)  -  1  for  all  e  >  0. 

3  2 

As  such,  from  (3.13),  we  get 

i  b— 

(3.16)  P(U(0  -  \f)  2  A-)  •  1. 

V°2 

Confining  (3.12)  and  (3.16),  we  get  P(U(U  -  W)*1  •  ~i --)  «  1 

Vb2 

which  in  turn  implies  U  »  bj,  V  ■  b^  Q 

Using  this  lemma,  we  can  establish  the  next  intermediate  theorem. 


1HB0BEM  3.3:  If_Yj  i  0,  and  if  are  independent,  then  VjYj  and  (1-V1)Y1  are 

dependent  if  and  only  if  one  of  the  five  conditions  below  are  true: 


SBBPf :  Th®  sufficiency  is  obvious,  so  we  need  only  establish  the  necessity  of  the 
conditions.  For  convenience,  let 

(3.15)  X  *  VjYj  ,  Z  -  (1-VjJYj. 

Further,  let 

Aj  «  {Yj  •  0}  o  {x  «  0,  Z  -  0) 

Aj  *  (Yj  >  0,  Vj  <  0}  »  {X  <  0,  Z  >  0} 

A3  »  {Yj  >  0,  Vx  «  0)  =  (X  ■  0,  Z  >  0} 

A4  «  (Yj  >  0,  0  <  Vj  <  1)  *  (X  >  0,  Z  >  0} 

Ag  ®  {Yj  >0,  Vj  s  1)  »  {X>0,  Z  •  0} 

A6  “  {Y1  >  °*  vj  >  1)  ■  (X  >  0,  Z  <  0). 

Also,  for  the  random  variables  Vj,  Yj,  X,  Z,  we  define,  respectively,  the 
variable  VA.  YA,  XA,  ZA  for  any  arbitrary  set  A  as  the  restriction  of  the  variable 
to  the  set  A.  That  is,  the  distribution  of  variable  VA  is  given  by 

(3.17)  P(VA  s  x)  -  [P(VeA)]”1  P(V  S  x,  VeA) 

with  the  other  variables  defined  identically. 

As  can  be  easily  established,  if  Vj  and  Yj  are  independent,  than  VA  and 

YAi  are  also  independent  for  all  i.  Similarly,  if  X  and  Z  are  independent  then  so 

are  XA.  and  Z,  for  all  i. 

Ai 

To  prove  the  theorem,  we  examine  three  cases, 

1)  P(Yj-0)  -  1  2)  P(Y1»0)  <  1,  P(0<VJ<1)  >  0  and  3)  P(Yj-O)  <  1,  P(0<V1<1)  -  0 

CASE  I:  P(Yj*0)  «  1, 

In  this  case,  we  have  nothing  to  prove,  this  being  condition  1)  of  the  theorem 


CASE  II:  P(Yj-O)  <  1,  P^Vj-cl)  >  0. 

In  this  case,  P(A4)  •  PCY^O)*  P^-eVj-el)  >  0. 

Clearly  X.  >  0,  Z.  >  0,  and  XA  and  Z.  are  independent.  Also  V.  ,  Y.  are 
A4  4  a4  a4  a4  4 

independent,  and  as  can  be  seen  from  (3.15),  we  have 

XA  4  ZA  *  YA  »  XA  (XA  ♦  YA  )_1  *  VA  * 

*4  n4  n4  *4  n4  "4  n4 

Thus  from  Lukacs  characterization  of  the  gamma  distribution  (see  [6]),  we  have  that 


(3.18) 


X.  -  r(A,a),  z.  ~  r(x,8) 

A4  a4 


Also,  we  have  that 


(3.19) 


YA  *  XA  4  ZA 
4  4  4 


Since  Y^y^}  ~  Y^  for  i  *  1  (from  the  independence  of  Yj  and  Vj),  we  have 
for  any  i  *  1  where  P(Ai)  >  0,  that 


(3.20) 


Y.  ~  Y.  for  i  a  1. 

Ai  a4 


By  similar  arguments,  we  can  also  establish  that  for  any  i  where  P(Ai)  >  0. 


we  must  have 


(3.21) 


X.  -  X.  for  i  >  4, 
Ai  A4 


Z  ~  Z  for  l<i<4. 
Ai  a4 


Let  us  assume  P^)  >  0.  Then  by  letting 


(3.22) 


0  *  \  W  * 


we  have  that  U  >  W  >  0.  Aseanbe  easily  established,  U,  df  are  independent,  and 
W(U-W)"1,  U-W  are  independent.  Also,  since  W(U-W)*1  -  U(U-W)"1  -l,  we  can  apply 


Lemma  3.2,  yielding  Z^  s  c,  X^  s  -d.  But  since  Z A  ~  T(A,B),  and  from  (3.21)  we 
have  Z  ~  ZA  ,  we  have  a  contradiction  which  yields  that  P(Aj)  «  0. 

If  P(Ag)  >  0,  we  can  generate  by  similar  methods  a  contradiction  for  the 
distribution  of  X.  ,  so  P(Afi)  ■  0.  If  P(A-)  >  0,  we  have  X.  ■  0,  so  Y,  ■  Z.  . 


rtX.o+6)  yielding  by  con¬ 


tort  from  (3.20)  snd  (3.21),  we  have  Y.  -  r(X,&)  and  Z,  ~ 

A3  a3 

tradiction  that  P(Aj)  *  0.  Similarly,  we  can  show  that  P(A^)  a  o,  yielding  that 
P(0<Vj<1)  •  1.  Thus  we  have  that,  letting  pj  ■  P(Yj  ■  0),  that 

(3.23)  P(XSx,  Zsz) 

-  Pl  I(X»0,Z*0)  ♦  (1-Pj)  T(X,a)Cx)  r(X,0)(z) 

Since  P(X*0)aP(Z«0)  «  P(X«0,  Z»0)  ■  ,  from  independence  of  X  and  Z,  we  have  that 

p^  ■  Pj2  which  yields  that  p^  «  0  since  p^  <  1.  This  yields  condition  5)  of  the 
theorem. 


CASE  III:  P(Y1=0)  <  1,  P(0<V.<1)  -  0. 

Since  P(Q<Vj<l)  «  0,  we  have  P(X>0,  Z>0)  ■  0.  Thus,  either  P(X>0)  «  0,  or 

P(Z>0)  «  0.  Assume  P(Z>0)  =  0.  Then  P(VjSO)  «  1.  In  addition,  assume  that 

P(V1<0)  >  0,  so  we  have  P(Aj)  >  0,  and  by  applying  Lemma  3.2  to  X.  ,  Z  ,  we  get 
Z^  *  c,  X^  •  -d,  c>d>0.  Thus  Y^  ■  Z^  ♦  X*  -  c-d>0.  If  PO^-O)  >  0  also, 

we  have  P(A-)  >  0.  By  (3.20)  and  (3.21),  we  get  that  Y.  ~  Y.  .  ?  c  -  d,  and 

s  a3  «2 

Z,  ~  Z.  »  c.  But  we  must  have  Z.  ■  Y,  .which  gives  a  contradiction.  Thus,  v 
.  a3  i  a3  o 

if  P(Vj  <  0)  >0,  then  P(Vj  •  0)  »  0  which  yields  that  P(V1  <  0)  ■  1. 


(3.24) 


Thus,  we  have  for  P(Yj  ■  0)  ■  Pj 
1  P(X-0,  z-0)  .  Pl 


[  P(X*-d,  Z-C)  -  (1-Pj)  , 


which  clearly  contradicts  the  independence  of  X  and  Z  unless  Pj  ■  0  or  *  1. 

Thus  we  have  X  and  Z  are  degenerate,  which  yields  condition  4). 

If  we  assume  instead  that  P(Vj»0)  >  0  in  addition  to  P(Z>0)  ■  0,  a  contra¬ 
diction  similar  to  the  one  above  will  show  that  P(Vj«0)  »  1,  which  is  condition  2). 


If  we  assume  that  P(X>0)  a  o,  arguments  identical  to  those  above  will  show 
that  P(Vj«l)  a  l  or  P(Vj>1)  =  1,  which  yield,  respectively,  condition  4)  or  con¬ 
dition  3).  □ 

We  are  now  ready  to  establish  the  main  theorem  of  this  section. 


THEOREM  3.4.  If  Vj ,  Yj  are  not  both  degenerate,  then  Vj ,  Y^^  independent  and  VjY^ 
and  (l-Vj)Yj  are  independent  if  and  only  if  one  of  the  five  conditions  below  are 
true: 

1)  Yj  5  0, 

2)  i  0, 

3)  Vj  s  1, 

4)  Yj  ~  r(X,04.p),  Vj  ** 

5)  -Yx  -  r(x,<**8).  Vj  ~  B(o,e). 


PROOF:  Again  the  sufficiency  is  obvious,  so  we  need  only  establish  the  necessity  of 


the  conditions.  If  we  have  either  PCY^O)  ■  1,  or  PCY^sO)  aj.t  then  by  theorem 


3.3  we  have  that  one  of  the  S  conditions  must  hold  (where  the  new  condition. 


condition  S)  follow  when  P(YjSO)  =1) .  Thus  from  now  on  we  will  assume 


PfY^O)  >  0,  PfY^O)  >  0. 


In 


addition  to  the  sets  to  Ag  defined  in  theorem  3.3,  we  define  the  sets 


A?  *  {Y^O,  Vj<0}  »  {X>0,  Z<0,  X+Z<0} 

A8  "  (Yj<°,  Vj-0)  *  {XaO,  Z<0) 

Ag  ■  (Y^O,  (>^<1}  -  {X<D,  Z<0> 

A10  *  {Y1<0»  Vl}  *  {X<0*  Z*0} 

^1  *  <Yi<0»  Vi>l)  8  {X<0»  Z>0»  x*z<0>» 


(3.25) 


and  we  redefine 


> 

^  »  {X<0,  Z>0,  X*Z>0} 
a6  (X>0,  Z<0,  X+Z>0} 


Note  that  while  for  all  i,  can  be  expressed  as  A^  ■  (YjtA,  V^eB)  (for  some 
borel  sets  A,  B)  the  sets  A^,  Afi,  Ay,  A^^,  cannot  be  expressed  as  sets  of  the 
fora  (XcA,  ZeB}. 

Let  us  assume  in  addition  to  P(Yj<0)  >  0,  P(Yj>0)  >  0,  that  P(0<Vj<1)  >  0 
(this  assumption  will  be  shown  to  yield  a  contradiction  later) . 

Thus  we  have  P(A^)  >  0,  P(Ag)  >  0.  By  appealing  to  Lukacs  characterization 
as  in  Theorem  3.3,  we  get 


(3.27) 


\  ~  r(Xl»al)*  \  ~  Ya4  ~ 

■XA9~  rtX2'a2?*  IP<W*  ‘YAg~  r<Xl»V0)‘ 


By  arguments  similar  to  those  used  for  (3.20)  and  (3.21),  we  have  that  if 
P(A6UA?)  >  0,  P(A2UAj1)  >  0,  then 


(3.28) 


\UA7*  ‘V'V"1'  *‘v''2’,'2' 

^AjUAjj  ^UAjj  ~  r(X2»®2)  x  r(Xl,0l)- 


~  r(A.,a1)  x  r(A,,B,) 


Since  P(AfiUA7)  >  0,  assume  that  P(A^)  >  0.  Thus  we  get  for  xQ  >  o,  yQ  <  0 


(3.29) 


P(XA  *  x0*  ZA  S  *Q>  ■  //— - 1 

*  0  A6  0  A  r(«1)r(B2)P(A6) 


°1  °l"^e”Xlx  02  X2Z 

X1  *  A2  *(-Z)  1  e  1  dxdz 


where  As((x,z)  :  x<xQ,  zSzQ,  x*z>0). 

Also,  by  arguments  similar  to  those  used  for  (3.20)  and  (3.21),  we  have  that 


-  19  - 


Thus,  by  transformation  of  variables  in  (3.28),  we  get  that 

6.  o.  ♦  6.-1  F(a..Bf)  1,62  -  a,-l  -U,*!,)*  8,-1 

<>•*>  h  y  1  •  Tisrj ikmj-to™  •  « 


By  examining  the  behavior  of  the  right  hand  side  of  (3.30)  and  the  left  hand  side 
of  (3.30)  (particularly  as  y  approaches  0),  we  get  that  equality  is  impossible 
for  >  0,  >  0.  Thus,  we  have  that  P(A^)  *  0.  Similarly,  we  can  establish 

that  P(A7)  ■  0.  Thus  P(X  >  0,  Z  <  0)  ■  P(A6UA?)  ■  0.  However,  since  we  have 
assumed  that  P(Yj  >0)  >0,  P(Yj  <  0)  >0,  and  P(0  <  Vj  <  1)  >0,  we  get  from  the 
definition  of  X  and  Z  (see  (3.15)),  that  P(X  >0)  >0,  P(Z  <  0)  >0.  This  yields 


a  contradiction  in  the  assumption  P(0  <  Vj  <  1)  >  0.  Thus,  we  conclude  that 
P(0  <  Vj  <  1)  ®  0,  and  we  have  that  P(X  >  0,  Z  >  0)  ■  P(X  <  0,  Z  <  0)  «  0. 

From  this,  we  can  quickly  deduce  the  rest  of  the  criteria.  Since 
P(Yj  >0)  >0,  and  P(Yj  <  0)  >0,  then  the  assumption  that  P(Vj  <  0)  yields  from 
(3.15)  that  P(X  -0)  >  0,  P(Z  >  0)  >  0.  Independence  of  X  and  Z  yields  that 
P(X  <  0,  Z  <  0)  >  0,  which  is  a  contradiction.  Thus,  P(Vj  <  0)  «  0.  In  a  similar 
manner,  P(Vj  >  1)  ■  0.  Thus  P(Vj  ■  0)  ♦  P(V^  »  1)  ■  1.  If  P(Vj  *  0)  >0, 

P(Vj  »  1)  >  0,  we  again  can  show  from  (3.15)  that  P(X  >0)  >0,  P(Z  >0)  >0, 
which  again  contradicts  P(X  >  0,  Z  >  0)  ■  0.  Thus  we  have  either  P(V^  «  0)  *0, 
P(Vj  ■  1)  »  1  (condition  3)  or  P(Vj  «  0)  el,  P(V^  ■  1)  ■  0  (condition  2).  Q 

To  relate  Theorem  3.4  to  the  asymptotic  independence  of  Yj,  n  and  Y2  n  , 
we  only  meed  to  observe  as  remarked  after  (3.7)  that  for  any  condition  sufficient 
to  guarantee  the  density  of  Ac  (see  (3.6)),  Yj  and  Y2  are  independent  if  and  only 
if  for  Vj,  Yj  independent,  V^Yj  and  (l-V^Yj  are  also  independent.  As  such,  for 
Ac  dense.  Theorem  3.4  gives  the  necessary  and  sufficient  condition  for 
asymptotic  independence.  The  most  reasonable  and  realistic  condition  on 
(Qn,Vn,Wn)  for  Theorem  3.4  to  be  applicable  are  given  in  the  following 
corollary. 


The  proof  follows  immediately  from  Theorem  3.4  and  the  fact  that  the 
characteristic  fractions  of  non-negative  random  variables  have  dense  support 


(see  Smith  (81). 

Clearly,  other  restrictions  to  C^,  VR,  and  will  yield  that  the  only  case 

of  asymptotic  independence  of  n  and  Y2  n  are  when  has  a  Ganna  distribution 

and  Vn  has  a  Beta  distribution.  However,  the  more  interesting  question  of  a 

characterization  of  the  distributions  of  0  ,  V  and  W  that  result  in  Y,  _  and 

ji  n  n  i#n 

Y_  being  asymptotically  independent  appears  to  be  an  open  question. 
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